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Abstract. We investigate the finite-time boundary stabilization of a 1-D first order quasilinear 
liyperbolic system of diagonal form on [0,1]. The dynamics of both boundary controls are 
governed by a finite-time stable ODE. The solutions of the closed-loop system issuing from 
small initial data in Lip([0,l]) are shown to exist for all times and to reach the null equilibrium 
state in finite time. When only one boundary feedback law is available, a finite-time stabilization 
is shown to occur roughly in a twice longer time. The above feedback strategy is then applied 
to the Saint- Venant system for the regulation of water fiows in a network of canals. 



1. Introduction 



Solutions of certain asymptotically stable ODE may reach the equilibrium state in finite time. 
iThis phenomenon, which is common when using feedback laws that are not Lipschitz continuous, 
was termed finite-time stability in [5] and investigated in that paper. 

A finite-time stabilizer is a feedback control for which the closed-loop system is finite-time 
stable around some equilibrium. In some sense, it satisfies a controllability objective with 
-a control in feedback form. On the other hand, a finite-time stabilizer may be seen as an 
J> exponential stabilizer yielding an arbitrarily large decay rate for the solutions to the closed- 
C^l loop system. This explains why some efforts were made in the last decade to construct finite- 
time stabilizers for controllable systems, including the linear ones. See |29[ [50] for some recent 
l/-^ developments and up-to-date references, and [2] for some connections with Lyapunov theory. 

• For PDEs, the relationship between exact controllability and rapid stabilization was inves- 
O .tigated in [311 1221 123] • (See also [24j for the rapid semiglobal stabilization of the Korteweg-de 
|Vries equation using a time- varying feedback law.) 

To the best knowledge of the authors, the analysis of the finite-time stabilization of PDE 
is not developed yet. However, the phenomenon of finite-time extinction exists naturally for 
certain nonlinear evolution equations (see [361 \TT\ [6]). On the other hand, it is well-known 
■since [28] that solutions of the wave equation on a bounded domain may disappear when using 
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'transparent" boundary conditions. For instance, the solution of the 1-D wave equation 



dh-dly = ^, in (0,r) X (0,1), (1.1) 

d^y{t,l) = -dty{t,l), in(0,r), (1.2) 

d^y{t,^) = dty{t,l), in(0,r), (1.3) 

(y(0,.),aty(0,.)) = (yo,^o), m(0,l), (1.4) 



is finite-time stable in {{y,z) G F^(0, 1) x ^^(0,1); 1/(0) + y(l) + /J z{x)dx = 0}, with T = 1 as 
extinction time (see e.g. |23[ Theorem 0.5] for the details.) The condition ()1.2p is transparent 
in the sense that a wave y{t,x) = f{x — t) traveling to the right satisfies ()1.2p and leaves the 
domain at x = 1 without generating any reflected wave. Note that we can replace (jl.Sp by 
the boundary condition y{t,0) = (or dxy{t,0) = 0). Then a finite-time extinction still occurs 
(despite the fact that waves bounce at a; = 0) with an extinction time T = 2. We refer to [8j for 
the analysis of the finite-time extinction property for a nonhomogeneous string with a viscous 
damping at one extremity, and to [1] for the investigation of the finite-time stabilization of a 
network of strings. 

The finite-time stability of (|l.ip - ()1.4p is easily established when writing (jl.ip as a first order 
hyperbolic system 

with (r, s) = {dxy, dty), and next introducing the Riemann invariants u = r — s,v = r + s that 
solve the system of two transport equations 

dtu + dxu = 0, 
dtv — dxV = 0. 

The boundary conditions (|1.2p and (|1.3p yield u{t, 0) = w(t, 1) = (and hence u{t, .) = v{t, .) = 
for t > 1), while the boundary conditions (|1.2p and y{t,0) = yield v{t,l) = and u{t,0) = 
v{t, 0) (and hence v{t, .) = for t > 1 and u{t, .) = for t > 2). 

The goal of this paper is to show that the finite-time extinction property can be realized for 
1-D first order quasilinear hyperbolic systems 

dtY + dxF{Y) = 0, (1.5) 

that can be put in diagonal form, i.e. for which there is a smooth change of (dependent) variables 
that transforms ()1.5p into a system of two nonlinear transport equations of the form 

dtu + A(u, v)dxU = 0, (1.6) 
dtv + ^{u,v)dxV = (1.7) 

where ^{u, v) < —c < c < X{u, v) are smooth functions and c > is some constant. In practice, 
the functions u and v are Riemann invariants of (jl.5p (see e.g. [13j). 

The generalization of the finite-time extinction property of the wave equation to systems of the 
form (jl.6p - (jl.7p is the main aim of this paper. Of course, one could just consider homogeneous 
Dirichlet conditions 

u{t,0) = v{t,l) = 0, 
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but this would impose to restrict ourselves to initial data (uo,fo) fulfilling the compatibility 
conditions 

uo(0) = voil) = 0. 

Rather, we shall consider boundary conditions whose dynamics obey a finite-time stable ODE, 
namely 

^n(t,0) = -Ksgn{u{t,0))\u{t,0)\\ (1.8) 

j^v{t,l) = -Ksgn{v{t,l))\v{t,l)P, (1.9) 

{{K,j) G (0, +00) X (0,1) being some constants) and supplement the system (|1.6|) - (|1.7|) . (|1.8p - 
with the initial condition 



u{0,x) = uq{x), v{0, x) = vo{x). (1-10) 

The first main result in this paper (Theorem[T]) asserts that for any pair (uq, vq) of (small enough) 
Lipschitz continuous initial data, system ()1.6p - ()1.7p and ()1.8p - ()1.10p admits a unique solution 
in some class of Lipschitz continuous functions, and that this solution is defined for all times 
t > and vanishes for roughly t > 1/c. Theorem [1] is proved by using a fixed-point argument 
(Schauder Theorem) and energy estimates. 

Sometimes, the boundary condition at one extremity of the domain (say 0) is imposed by the 
context, so that we cannot chose the condition n(t,0) = (or its generalization (|1.8p ) for the 
Riemann invariant u. Then, we have to replace (|1.8p by a boundary condition of the form 

u{t,0) = h{v{t,0),t), (1.11) 

for some (smooth) function h = h{v,t). The second main result in this paper (Theorem [2]) 
asserts that the system ()1.6p - (|1.7p and ()1.9p - ()l.lip is still locally well-posed with roughly an 
extinction time T = 2/c. The result is obtained for small initial data and for ||c?(/i||oo small 
enough. 

The results obtained in this paper can be applied to: 

(1) the p— system 

dtr - d^s = 0, (1.12) 
9iS-9,[p(r)] =0, (1.13) 

where p G C^(M) is any given function; 

(2) the shallow water equations (also called Saint- Venant equations [55] ) 

dtH + d^{HV) = Q, (1.14) 

dtV + d^A— + 9H) = Q, (1.15) 

where H is the water depth and V{t^ x) the averaged horizontal velocity of water in a 
canal, and g the gravitation constant; 

(3) Euler's equations for barotropic compressible gas 

dtp + d^{pV) = Q, (1.16) 
dt{pV) + d^{pV^+p) = {), (1.17) 
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where p is the mass density, V the velocity, and p = p{p) the pressure of the gas. 
(4) The same strategy could in theory be applied to any system possessing Riemann invari- 
ants. Riemann invariants exist for most 2x2 systems, and also for some larger systems 
(e.g. the 3x3 system of Euler's equations for compressible gas, see [35, chapters 18,20]). 

For the sake of shortness, we will limit ourselves to the stabilization of Saint- Venant equations, 
and will give an extension of the above finite-time stabilization results to a tree-shaped network 
of canals. The obtained extinction time will be roughly d/c, where d denotes the depth of the 
tree (Theorem [5]) . 

There is a huge literature about the controllability and stabilization of first order hyperbolic 
equations (see e.g. [151 El ESI EH [271 EHl HH]). In particular, the control of Saint- Venant 
equations has attracted the attention of the control community because of its relevance to the 
regulation of water flows in networks of canals or rivers. We refer the reader to e.g. [3 [571 [^ 
[2CT1 [171 [11 [T?l [1^1 [3], where Riemann invariants played often a great role in the design of the 
controls. Our main contribution here is to notice that a finite-time stabilization can be achieved 
as well, i.e. that bounces of waves at the two ends of the domain can be avoided. 

A numerical scheme and some numerical experiments for the finite-time stabilization of water 
flows in a canal may be found in [32], in which certain results of this paper were announced. 

The paper is outlined as follows. Classical but important properties of linear transport equa- 
tions are recalled in Section 2. In Section 3, we introduce two boundary controls whose dynamics 
are governed by a finite-time stable ODE, and prove the existence and uniqueness of a solution 
to the closed-loop system, and the fact that this solution reaches the null state in finite time. In 
Section 4, we investigate the same problem with only one boundary control, the other boundary 
condition being imposed by the physical context. In the last section, we apply the results in 
Sections 3 and 4 to the regulation of water flows in a canal with one or two boundary controls, 
and extend the finite-time stabilization results to any tree-shaped network of canals. 

2. Some background about linear transport equations 

2.1. Notations. e°([0,r] x [0, 1]) denotes the space of continuous functions u : [0,T] x [0, 1] — ^ 
M. It is endowed with the norm 

ll'"lleo([o,T]x[o,i]) = sup \u{t,x)\. 

{t,x)e[o,T]x[o,i] 

The norm of the space LP{0, 1) is denoted || • ||p for 1 < p < oo. Lip([0, 1]) denotes the space 
of Lipschitz continuous functions u : [0, 1] — )• M. It may be identified with the Sobolev space 
Wi'°°(0, 1). Lip([0, 1]) is endowed with the W^'°°{0, l)-norm; that is 

ll^llLip([0,l]) = ll'"llvi/i>°°(0,l) = Halloo + IK'Hoo- 

We use similar norms for Lip(M), Lip([0,T] x [0, 1]), etc. 

2.2. Linear transport equation. In this section we consider the initial boundary-value prob- 
lem for the following linear transport equation 



dty + a{t,x)dxy = 0. 



(2.1) 
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We assume thereafter that 

a e e°([0, T] X [0, 1]) n L°°(0, T; Lip([0, 1])), (2.2) 

a{t,x)>c>0, V(t,x) G [0,r] X [0, 1], (2.3) 

where c denotes some constant. Note that the case when a{t,x) < —c < can be reduced to 
()2.3p by the transformation x — t- 1 — x. 

2.3. Properties of the flow. By (|2.2p . a is uniformly Lipschitz continuous in the variable x, 
with say a Lipschitz constant L = ||a||/^oo(o,T;Lip([o,i]))- Since we intend to use the method of 
characteristics to solve (|2.1|) . we need to study the flow associated with a. 

Definition 1. For {t,x) G [0,T] x [0,1], let (/){., t,x) denote the maximal solution to the 
Cauchy problem 

ds(l){s,t,x) = a{s,(l){s,t,x)), ^^^^ 
(p{t, t, x) = X, 

which is defined on a certain subinterval [e{t,x), f{t,x)] of [0,T] (which is closed since [0,1] is 
compact), and with possibly e{t,x) and/or f(t,x) = t. Let 

Dom(l) = {{s,t,xy, (t,a;) G [0,T] X [0,1], s G [e{t,x), f{t,x)]} 

denote the domain of (p. 

Note that 

e{t, x) > ^ (j){e{t, x),t, x) = 0. (2.5) 
We take into account the influence of the boundaries by introducing the sets 
P := {{s,^{s, 0,0)); s G [0, /(O, 0)]}, 
/ := {(t,x) G [0,T] X [0,1]\P; e(t,x) =0}, 
J := {(t,x) G [0,r] X [0,1]\P; ,/.(e(t,x),t,x) =0}. 
(See Figure [U) Note that both / and J are open in [0,T] x [0, 1]. 




e{t.x) T t 



Figure 1. Partition of [0,T] x [0, 1] into / U P U J. 
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Proposition 2.1. Let a satisfying ()2.2p . let L = | |a| |Loo{o,T;Lip{[o,i])); '^^^ 

K := max(l, ||a||eO([o^T.]x[o,i]))e^'^- 

Then (f) is K-Lipschitz on its domain; that is, for all (si,ti,xi), {s2,t2,X2) G Domcf) 

|(/)(si,ti,Xi) - (t){s2,t2,X2)\ < K{\si - S2I + \tl - t2\ + \xi - X2I) • (2.6) 

The proof of Proposition 12 . 1 1 is given in appendix, for the sake of completeness. 
We can now study the regularity of e. 

Proposition 2.2. Let a be as in Proposition\2^ let {t,x) G [0,T] x [0, 1], let {un} C e°([0,r] x 
[0, 1]) n L°°(0, T; Lip([0, 1])) be a sequence such that | |an| |L°o(o,T;Lip([o,i])) ^■s bounded and 

\\an - a||eO{[o,T]x[o,i]) -^0 as +oo, 
and let {{tn,Xn)} C [0, T] x [0, 1] be a sequence such that {tn,Xn) {t,x). Then 

en(tn,Xn) ^ e{t,x). (2.7) 

Proof. We use again the extension operator 11 introduced in the proof of Proposition 12.11 (see 
the appendix) and set a„ = n(a„) and d = n(a). Let (pn and (p denote their respective flows. 
Recall that (j) and (j) coincide on Dome/) (resp. and (pn coincide on Dom (/;„). From 

\ds[<i)n{s,t,x) - (i){s,t,x)]\ = \dnis,(f)n{s,t,x)) - d{s,(j){s,t,x))\ 

< \dnis,(t>n{s,t,x)) - d{s,(t>nis,t,x))\ 

+ \d{s,(j)n{s,t,x)) - d{s,(f>{s,t,x))\ 

< ll^n - a||Lcx>(K2) + ||a||ioo(iR.Lip(iR))|(^„(s,t,x) - 4>{s,t,x)\, 

(fOnj) and Gronwall's lemma, we infer that for all n > and all is,t,x) G [0,T]^ x [0, 1], 
we have 

\{4>n - 4>){s,t,x)\ < T\\an - a||eo([o,T]x[o,i])e''ll'^ll^°°(°™(M». (2.8) 
It may be seen that 

e„(i„,x„) = min{s G [0,t„,]; VrG[s,t„], (^„(r, t„, x„) G [0, 1]}. 

• If (t,x) G /, then since we have excluded the characteristic coming from (0,0), we have 
that 

inf dist((s,(/)(s,i,x)), [0,t] x {0}) > 0, 

sG[0,i] 

where dist((t, x), F) = inf(j/ .j,/)gp'(|t — t'\ + \x — x'\). So we infer from ()2.8p that for n 
large enough (j)n{-,t,x) is defined on [0,t], i.e. en(t, x) = 0. Then (|2.7|) is obvious. 

• From now on, we assume that (i, x) G J L) P. We claim that 

limsupe„(t„,x„) < e{t,x). (2.9) 

Indeed, if e{t, x) = t, then 

limsupe„(t„,Xn) < limsupt„ = t = e{t,x). 

n— >oo n—>ca 



FINITE- TIME STABILIZATION OF SYSTEMS OF CONSERVATION LAWS ON NETWORKS 7 

Otherwise, we have e{x,t) < t, and using ()2.3p we obtain for any e E (0, (t — e{t,x))/2), 
ce < (j){s,t,x) <1 - ce, Vs E [e(t, x) + e, t - e]. (2.10) 
However, we have for n large enough 

~ ~ C€ 

\\4>n - </'lleO([0,T]2x[0,l]) < J> (2-11) 

~ ~ C6 

\(pn{s, tn, Xn) - (t>n{s, t,x)\ < — , (2.12) 

the second estimate coming from the uniform bound on ||an||L°o(o,T;Lip([o,i])) Propo- 
sition 12.11 Combining ()2.10p . ()2.1ip and ()2.12p . we see that for n large and for ah 
s E [e{t, x) + €,t — e], (pnis, tn, Xn) is well defined and 

ce 

(PnySjtm Xn) ^ ~- 

This yields lira supn^oo ^n{tm ^n) < e{t,x) + e, and since e was arbitrarily small, (|2.9p 
follows. 

If {t,x) E P the proof of (j2.7p is complete, for lim inf„^oo en(tn, 2;„) > = e{t,x). 
Assume finally that (t, x) E J, so that e{t, x) > 0. Pick any s E (0, e(t, x)). We obviously 
have (p{s, t, x) < 0, thanks to the lower bound on o (see ()5.4ip ). But we know from ()2.8p 
and Proposition 12.11 that 

4>n{s,tn,Xn) 4>{s,t,x), 
n^-\-oo 

and hence for n large enough, (pn{s,tn, Xn) < and s < e(tn,Xn)- Thus, we conclude 
that liminf„_!.oo en{tn,Xn) > s. As s was arbitrarily close to e{t,x), we end up with 

liminf e„(t„,x„) > e{t,x). 

n— >oo 

The proof of (j2.7p is complete. □ 

Remark 1. (1) For Un = a, this shows that e is continuous on [0,T] x [0, 1]. 

(2) Since [0,T] x [0, 1] is compact, Proposition \2.^ implies that Cn converges uniformly toward 
e on [0,T] x [0,1]. 

Proposition 2.3. //, in addition to <^2!l-<^M' we have that d^a E e°([0,r] x [0, 1]), then (j) is 
on Dom (p and e is on [0, T] x [0, 1] \ P, with for {t, x) E J 

a{t,x)exp{- f\^Sa:a{r,(t){r,t,x))dr) exp{~ j\^Sxa{r,(f){r,t,x))dr) 

dte{t,x) = \ , d.j,e{t,x) = \ 

a(e(t,x),0) a(e(t,x),0) 

(2.13) 

Proof. The regularity of (/> is a classical result (see e.g. [2T]). If (t,x) E /, e(t,x) = and the 
result is obvious. For {t, x) E Jn (0, T) x (0, 1) we have (l){e{t, x), t, x) =0 and ds(j){e{t, x),t, x) > 
0, therefore the Implicit Function Theorem allows us conclude. Finally, for (t, x) E J \ (0,T) x 
(0, 1), it is sufficient to pass to the limit in ()2.13p . □ 
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Proposition 2.4. Let a fulfill ()2.2p and ()2.3p . and let L = ||a||Loo(o,T;Ljp([o,i]))- Then the 
function e is K -Lipschitz on [0, T] x [0, 1] where K is given by 

K = c"^ max {l, ||a||eo([o,T]x[o,i]))e'^^- 

Proof. Consider (ti, xi) and {t2^X2) in [0, T] x [0, 1]. Let us also suppose tliat e(ti, xi) > e(t2, 3^2), 
tlie other case being symmetrical. We infer from Proposition 12.11 that 

|(/)(e(ti,xi),ti,xi) - (j){e{ti,xi),t2,X2)\ < max (^1, ||a||eo([o,T]x[o,i]))e^^(|ii - hi + \xi - X2I). 

(2.14) 

Since e(ti,xi) > 0, we have that (j){e{ti,xi),ti,xi) = 0, and 

(pie{ti,xi),t2,X2) > c(e(ti, xi) - e{t2,X2)) > 0. (2.15) 
Therefore we end up with 

|e(ti,xi) -e(t2,X2)| <K{\ti-t2\ + \xi-X2\). (2.16) 

□ 

2.4. Strong solutions. Let a £ e°([0, T]; e^([0, 1])), yi G e^([0,r]), and yo G £^[0,1]) be 
given, and assume that the following compatibility conditions hold: 

ym = yo(0), yRO) + a(0, 0)y^,(0) = 0. (2.17) 

We consider the following boundary initial value problem: 

dty + a{t,x)d^y = 0, (t, x) G (0, T) x (0, 1), (2.18) 
y{t,0) = yi{t), tG(0,r), (2.19) 
y{0,x)=yoix), xG(0,1). (2.20) 

A strong solution of (IXTB|) - (p:^ is any function y G e'^{[0,T] x [0,1]) such that (IXTB|) - (IX^ 
hold pointwise. 

We define a function y : [0, T] x [0, 1] — )• M in the following way: 

^^ = if(i,x)GJ, 

' ' \yo{<PiO,t,x)) if (t,x) GIUP. ^ ■ ^ 

Proposition 2.5. Let y 6e as in ()2.2ip . T/ien y is a strong solution of ()2.18p - ()2.20p . Besides, 
we have the estimates 

l|y|leO([o,T]x[o,i]) < ™ax (||yo||eo([o,i])) llydleO([o,T])) > (2.22) 

IIV7 II ^ /^ll 'II II 'II ll«lleO([0,T]x[0,l])\ II N 

II vy||eo([o,T]x[o,i]) < max I ||yolleo([o,i])) llyilleo([o,r]) I exp [I \\Oxa\\Qo^io^T]x[o,i])) ■ 

(2.23) 
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Proof. One can see that y is of class on I and J, with the derivatives given by: 

dty{t, x) = y'i{e{t, x)) exp ( - / da:a{s,(l){s,t,x))ds] , V(t,x) G J, 

a{e{t,x),0) y J^^t,x) J 

dxy{t,x) = -y[{e{t,x))—-—^———exp i- [ dxa{s,(j){s,t,x))ds] , V(t,x) G J, 

a{e{t,x),0) y J^^t^^) J 

dty{t, x) = -2/0(0(0, t, x))a{t, x) exp ^- j d^ais, (j){s, t, x))ds^ , \/{t, x) G /, 

d:ry{t,x) = yQ{(f>{0,t,x))exp j d^a{s,<f>{s,t,x))ds^ , \/{t,x)eI. 

It follows from the first equation in ()2.17p and the continuity of e that y is continuous at each 
point of P. Note that y is differentiable in directions t and x in the following way: for all 
i G (0,/(0,0)) 

1 =yi — — exp - / 9^a(s,(/.(s,0,0 )ds , 

/i^o+ Ai a(0, Oj \ Jo J 



lim ^(^ + ^'^(^'0'")) = _y^(o)a(t,0(i,O,O))exp f- [' dMs,^is,0,0))ds 

hm = 2/0(0) exp - / d^a{s, (j){s, 0, 0))ds , 

h->o+ n. V Jo / 

hm = -y;(0 — — exp - / 9^a(s, (/>(s, 0, jds 

,^0- ft a(0,0) V Jo 



h^o- 

Using the second equation in ()2.17p , we see that y G Q^{[Q,T] x [0, 1]). The fact that y satisfies 
()2.ip follows from a straightforward calculation. □ 

2.5. Weak solutions. Now we consider the case when a G e°([0,r] x [0,1]) and d^a G 
L°°((0,r) X (0, 1)). We stih assume that 

a(i,a;)>c>0, V(t, x) G [0, T] x [0, 1]. (2.24) 

We begin by introducing the space: 

T={V'Gei([0,r] X [0,1]); V(i,l) = V'(r,x) =0 V(t,x) G [0,T] X [0,1]}. (2.25) 

We say that a function y G -^^^(0) T) x (0, 1)) is a weak solution of (|XTS|) - ([X^ if for any V' G T 
we have 

y(t, x){tpt{t, x) + a(t, x)ipx{t, x) + aa;(t, x)ip{t, x))dtdx 

(0,T)x(0,l) 

+ / V(i,0)y/(t)a(t,0)(it + 7p{0,x)yo{x)dx = 0. (2.26) 



Using the results of Section 12.41 it is clear that a strong solution is also a weak solution. Con- 
versely, any weak solution which is in C"'^([0,T] x [0,1]) is a strong solution. Note that the 
definition of weak solution makes sense for yi G L^{0,T) and yo G 1^(0, 1). 
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Proposition 2.6. Let us suppose that a, yi and are uniformly Lipschitz continuous with 
Lipschitz constants L, Li and Lq, respectively, and that yi{0) = yo(0)- Then the function y 
defined by 

' ' \yoimt,x)) if{t,x)eIUP, ^ ' ' 

is a weak solution of (j2.18p - (j2.20p . Furthermore, y is M-Lipschitz continuous on [0, T] x [0,1] 
with M defined by 

M := max(^,Lo)max (l, ||a||eo([o,T]x[o,i]))e'^^- (2-28) 

Finally, y is the unique solution in the class Lip{[0,T] x [0,1]) of system (|2.18p - (|2.20p . with 
(|2.18p understood in the distributional sense, and (j2.19p - (|2.20p pointwise. 

Proof. Using standard regularization arguments, it is possible to find a" G C^([0,T] x [0,1]), 
G ei([0, 1]) and G ei([0, 1]) such that: 



Ve E (0,min(r, 1)) ||a'' - a||eO([o,T]x[o,i]) + l|yo - yolleO([e,i]) + ll^r - yi\\&{[e,T]) 0, 

(2.29) 

ll«n||eO([0,T]x[0,l]) < 2||a||eO([o,T]x[0,l])i l|f^xOn||L°°((0,T)x(0,l)) < 2| |9a;0| |ioo((o,T)x (0,1)) > (2.30) 

l|yolleo([o,i]) < 2||?/o||eo([o,i])> ll?/o'lk°=(o,i) < 2||yollL°°(o,i), (2.31) 
l|yrileo([o,T]) < 2||yi||eo([o,T])) l|yr'llL°°(o,T) < 2||?/;||loo(o,t), (2.32) 

yo(o) = yr(o), yo"'(o) = yr'(o) = VneN. (2.33) 

Using Proposition 12.51 we infer the existence of a strong solution G C^([0,T] x [0, 1]) of 

fa,y" + a"9.y« = 0, V(t, x) € (0, T) x (0, 1). (2.34) 

\y-(i,0)=yr(t), y^{Q,x) = y-{x), 

y"- is given by ()2.2ip . with yi,yQ,e and (p replaced by y",yo'^'^ ^"^-i respectively. Note that 
(t, x) € (resp. (t, x) G J") for n large enough if {t, x) £ I (resp. {t, x) £ J). Using Proposition 
EH dSS]) and ^TTDf . we see that 

y"(t,j;) y{t,x), V(i,x)eIUJ. (2.35) 

n— >+oo 

Note that P = [0, T] x [0, 1] \ (/U J) has zero Lebesgue measure. An application of the dominated 
convergence theorem yields 

-y|lLi((o,T)x(o,i)) 

Using the other convergence assumptions about o", yf, and t/q, we can pass to the limit in 
(|2.26p . This shows that y is a weak solution of (|2.ip . 

To prove the regularity of y we distinguish two cases. 
Assume first that both (ti,xi) and {t2,X2) are in J U P. Using ()2.2ip and Proposition 12.41 we 
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have that 

\y{ti,xi) - y{t2,X2)\ = |yz(e(ti,a;i)) - yi{e{t2,X2))\ 
< Li\e{ti,xi) - e{t2,X2))\ 
L 



< 

c 



-max(^l,||a||eo([o,r]x[o,i]))e^^(|ii -^2! + \xi - X2\). 



Next, if we assume that {ti,xi) and {t2, X2) are in /U P, then we can use (|2.2ip and Proposition 
l2.1l to obtain that 

\y{ti,xi) - y{t2,X2)\ = \yo{4>{0,ti,xi)) - yo{4>{0,t2,X2))\ 

< Lo|0(O,tl,Xi) -0(O,t2,X2)| 

< Lo max (l, ||a||eo([o,T]x[o,i]))«^^'^(l*i " ^2! + \xi - X2I). 

Finally, if {ti,xi) E J and (t2,X2) G /, we consider an intermediate point on P belonging to 
the boundary of the rectangle [mm(ti,t2),max{ti,t2)] x [min(3;i, 3:2), max(2;i, X2)] and use the 
estimates above. 

Let us now check that y is the only solution to (|2.18|) - (|2.2U|) in the class Lip([0,T] x [0, 1]). 
First, picking any ijj G eg°((0, T) x (0, 1)) in I^TM . we see that ([21^ holds in D'((0, T) x (0, 1)). 
Note that each term in belongs to L°°((0,T) x (0, 1)), so that ([XTH|) holds also pointwise 

a.e. Scahng in ([TTH|) by il^ e 7 and comparing to (fTM]) . we obtain that (^1^ and ([X^ 
hold a.e., and also everywhere by continuity of yi, yo, and y. Thus y solves (|2.18|) - (j2.2U|) . If 
y G Lip([0,r] X [0, 1]) is another solution of (IXTH|) - (|CT]) . then y ■= y - y e Lip([0,T] x [0, 1]) 
solves 

dty + a{t,x)d^y = in a)'((0, T) x (0, 1)), (2.36) 
y{t,0) = in(0,T), (2.37) 
y(0,x) = in (0,1). (2.38) 

Scaling in ()2.36p by 2y, integrating by parts and using ()2.3p . ()2.37p . and ()2.38p . we obtain 

\\m\\l= f!\d^a)\y\^dxds- f a{t,l)\y{t,l)\^dt < L f \\y{s)\\lds. 
JoJo Jo Jo 

This yields y = by Gronwall's lemma. The proof of Proposition 12.61 is complete. □ 

3. Finite-time boundary stabilization of a system of two conservation laws 
In this section, we consider the system 

* . ' (t,x)e(0,+oo)x(0,l). (3.1) 

dtv + n{u,v)dxV = 0, 

where A and /i are given functions with 

A,/i G e°°(R^M), (3.2) 
/i(u, w) < -c < < c < X{u, v), y{u, v) G (3.3) 
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for some constant c > 0. We aim to prescribe a control in a feedback form on the boundary 
conditions u{t, 0) and v{t, 1) so that for some time T we have for any small (in Lip([0, 1])) initial 
data uq and vq 

u{T,x) = v{T,x) = 0, VxG(0, 1). (3.4) 

Remark 2. (1) If we intend to stabilize the system around a non null (but constant) equi- 
librium state {u,v) G M?, it is sufficient to consider the new unknowns u := u — u, 
V := V — V that satisfy a system similar to ()3.ip . and to stabilize {u,v) around (0,0). 
(2) Note that, since we are only interested in proving a local stabilization result, the condition 
()3.3p is not too much restrictive. It should be seen as X{u,v) > and fj,{u,v) < 0. 



After introducing the boundary feedback law, we will show the existence and uniqueness of 
the solution to the closed loop system and check that the property ()3.4p indeed holds for this 
choice of feedback law. 

We now come back to the quasilinear system (j3.ip that we complete as follows: 



dtu + A(u, v)dxU = 0, (t, x) G (0, +oo) x (0, 1), 
dtv + fj.{u, v)dxV = 0, (t, x) G (0, +oo) x (0, 1), 

^u{t,d) = -Ksgn{u{t,m<tM\ t > 0, 
at 

^^v{t,l) = -Ksgn{v{t,l))\v{t,l)\\ t > 0, 



(3.5) 



(3.6) 



n(0, x) = no(x), v{0,x) = vo{x), x G (0, 1) (3.7) 

with (-fC, 7) G (0,+oo) X (0,1) arbitrarily chosen. We aim to use Schauder fixed-point theorem 
to prove the local in time existence of solutions {u,v) of (j3.5p - (|3.7p in some class of Lipschitz 
continuous functions. By solution, we mean that (j3.5p is satisfied in the distributional sense, 
and that (|3.6p - (|3.7p are satisfied pointwise. Actually, we shall use the results of the previous 
section and define u as the weak solution of the transport equation (j'2.18p - (|'2.20p with a{t,x) = 
\{u{t,x),v{t,x)) for some given pair {u,v) in the same class, yi{t) = ui{t) (see below (|3.10p ). 
and yo{x) = uq{x), and similarly for v. 

3.1. Notations. Let Ci > and C2 > be given, and pick any uq,vq G Lip([0, 1]) with 

max(||-uo||oo, Ibolloo) < Ci, (3.8) 
max(||n[)||oo,|bol|oo) <C2. (3.9) 

Let 

r:=-+ ^1 



c {l-l)K 

We define ui and as the solutions of the following ODEs 

d 



uiit) = -Ksg-a{ui{t))\ui{t))W ui{Q) = no(0), 

(3.10) 

-gVr{t) = -KsgIi{Vr{t))\Vr{t)\\ Vr{^) = Vq{1). 
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An obvious calculation gives 

\uo{0)l~-' 



Ul{t) = < 

and 

Vr{t) = < 

Clearly 



sgn(no(0)) (|no(0)|i-^ - {I - ^)Kt)— if < t < 
ifi>M0^- 



(i-7)i<: ' 



sgn(^o(l)) (bo(l)l^-^ - (1 - l)Kt)^< if < t < ^l^iill^, 

"o(i)r-^ 

{l-y)K ' 



ift>MiriZ 



Vt>r--, vM) =ui(t) = i), (3.11) 
c 

max(||ni||oo, ll^^rlloo) < C*!, (3.12) 

max(||uJ||oo, ||fr||oo) < -^^(^7. (3.13) 

Let us also introduce 

Ml := max (||A||eo([_Ci,Ci]2). ll/^lleo([-Ci,Ci]2)) > (3-14) 
M2 := max (^\^u^l\\&o^^^_c^^Cl]■^)^\^y^l\\QOl^\^_c^^Cl]^^)^\^uA^^ 

(3.15) 

Let us pick a positive number C3. Let T> denote the domain 

D := |(ti,w) € Lip([0,T] x [0,1])^; max(||n||oo, ||i^||oo) < C*!, and u and w are Cs-Lipschitzj. 

(3.16) 

Let us equip the domain T> with the topology of the uniform convergence. Then, by Ascoli-Arzela 
theorem, 2) is a compact set in e°([0,r] x [0,1])^. 
The main result in this section is the following 

Theorem 1. Assume that Ci > and C2 > are such that 

/ \ /KC^ \ 1 
TM2 max 1, Ml max ^, C2 < — (3.17) 



c J - 2e 

and let C3 = {2TM2)~^. Pick any pair (uq^vq) £ Lip([0, 1])^ satisfying (j3.8p - (j3.9p . Then there 
exists a unique solution {u,v) of ()3.5p - ()3.7p in the class D. Furthermore, the solution is global 
in time with u{t,.) = v{t,.) = for t > T. Finally, the equilibrium state (0,0) is stable in 
Lip{[0,l]f for that is 

ll(«,''^)llL-'(R+;Lip([0,l])2) ^ as ||('U0,1'0)||Lip([0,l])2 ^ 0. (3.18) 

The first task consists in constructing a solution of the closed loop system as a fixed point of 
a certain operator. 
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3.2. Definition of the operator. If (tt, v) G T) are given, we define (n, v) = 3'{u, v) as follows: 
the function u is the weak solution of the system 



{dtu + A(n, v)dxU = 0, 
u{t,0) = ui{t), u{0,x) = uo{x), 

and the function v is the weak solution of the system 

idtv + iJ.{u,v)dxV = 0, 

\v{t,l) = Vr{t), v{0,x) = Vo{x), 



V(t,x) G [0,r] X [0,1], (3.19) 



V(t,x) G [0,r] X [0,1]. (3.20) 



3.3. Stability of the domain. In this part, we show that for a certain choice of Ci, C2, C3, 

we have 

^(03) C V. 

We first apply the results of Section 2 to get the following 

Lemma 1. Let Ci, (72,6*3 be any positive numbers, and let uo,vo G Lip{[0,l]) satisfying (|3.8p - 
()3.9p . For ^'ii'en (n,-?) G D, /ei (n, t>) = 3'{u,v). Then the functions u and v are Lipschitz 
continuous on [0, T] x [0, 1] and they satisfy the following estimates 

max(||n||eo([o,T]x[o,i]) JblleociCTlxfo,!])) < Ci, (3-21) 

max(||92:n||ioo((o^T)x(Oa))' ll^2^^lli°°((0,T)x(0,l)) J|t^i^llL°°{(0,T)x(0,l)) J|t^t^llL°°({0,T)x(0,l))) 

< max ^, C2) max (l, Mi) exp (^2TM2C3) . 

(3.22) 

Proof Estimate l^7II\i follows directly from (|CTD . (|CTD . (g^ZD, (I2S1) and (fXT^ . 

Estimate (IH:^^ can be deduced applying (IX^ for (IXTUD and (IX^Uj) . and using ([MD, ([HT^ . 
(fmil . (fXT^ and dSISD- □ 

Thanks to Lemma [TJ we see that the domain D is stable by 3" as soon as 
max f 



(^^,62) max (1, Ml) exp (2rM2C3) < C3. (3.23) 



This can be written as 



/KC^ \ C'3exp -2rM2C3 

max ( C2 ) < ^ r ^. (3.24) 

^ c / max (1, Ml 



For given Ci and C2, T, Mi and M2 are fixed. Note that T, Mi and M2 are independent of 
C2, and that they are nondecreasing in Ci. Therefore, as a function of C3 the supremum of the 
right-hand side of ()3.24p is attained for C3 = (2TM2)~^, and for this value of C3 the condition 
on Ci and C2 for the domain to be stable reads 

/ \ /KC^ \ 1 
rM2 max I 1, Ml ) max ( ^ , C2 I < — . (3.25) 



c / - 2e 

But the term in the left-hand side of (j3.25p tends to when Ci and C2 tend to 0, so that for 
Ci , C2 small enough the condition ()3.25p is satisfied and D is stable by 3". 



FINITE- TIME STABILIZATION OF SYSTEMS OF CONSERVATION LAWS ON NETWORKS 15 

3.4. Continuity of the operator. In this part we consider a sequence {{un,Vn)} C 2) and a 
couple {u, ?}) G D such that 

max (^\\un - u\\eo([o,T]x[o,i])^ W^n - ^'lleo([o,r]x[o,i])) ^J^^ 0. (3.26) 

Let us now define 

{un,Vn) = 3^{un,Vn) for n > 0, and {u,v) = 3'{u,v). (3.27) 
Our goal in this subsection is to show that 

max (^\\un - u\\qo(^[o^t]x[o,i])^ W^n - f^lleo([o,T]x[o,i])) 0- (3.28) 

We need the following 

Lemma 2. For almost all {t,x) S [0, T] x [0,1], we have 

{Un{t,x),Vn{t,x)) {u{t, x) ,v{t, x)) . (3.29) 

n— ^+00 

Proof. Let us show that Un{t,x) — t- u{t,x), the convergence Vn{t,x) — )• v{t,x) being similar. 

The fact that (u„, {;„) converges uniformly toward {u, v) on [0, T] x [0, 1] implies that \{un, Vn) 
converges uniformly toward X{u,v) on [0,T] x [0,1]. Furthermore, since (?!„,?}„) G T) for all n, 
we see that the functions A('u^, Vn) are uniformly Lipschitz continuous for n > 0. This will allow 
us to use Proposition 12.21 To this end, we consider the flow cpn (resp. cp) of A(n„,t}„) (resp. 
X{u,v)). In the same way, we define e„ and e, I„ and /, J„ and J, P„ and P. Using ()2.27p we 
have that 

nn{t,x) = h'-^'/''^\^ if(t,-)eJ„, 

\uo{cl)n{0,t,x)) if{t,x)eInUPn, 

and also 

^)^lMe{t,x)) ii{t,x)eJ, 
^' ^ |uo((/>(0,t,x)) if(t,x)G/UP. ^' ^ 

We infer from Proposition 12.21 that 

en{t,x) e{t,x), V(t,x) G [0,T] X [0, 1]. (3.32) 

This shows in particular that if {t,x) G J, then e{t,x) > and hence en{t,x) > for n large 
enough, i.e. (i, x) G Jn for n large enough. Therefore 

Un{t,x) — )• u{t,x), V(t, x)GJ. 

n— >-(-oo 

Now if (t, x) G /, then e(t, x) = and 0(0, t, x) > 0. Since A > c > 0, this implies the existence 
of e > such that 

e < (l){s,t,x), VsG[0,t]. (3.33) 

Combined with (|2.8|) . this shows that for n large enough e„(t, x) = and (/>ri(0, t, x) — )• (p{0, t, x), 
so we conclude that 



Un{t,x) = UQ{(j)n{0,t,x)) -> Uo{4i{0, t, x)) = u{t, x) . (3.34) 

71— > + 00 

Finally, P is clearly neghgible and / U P U L = [0, T] x [0, 1]. □ 



max 
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To strengthen this convergence, we just need to recah that for every n > 0, we have v„) E 
D which is compact in C'^([0, T] x [0,1]). According to Lemma [21 the only possible limit point 
is {u,v) and therefore we get the convergence of the whole sequence in 2); that is, 

(\\Un - u\\eO{[0,T]x[0,l])A\Vn - ^^1 leO([0,T] X [0,1])) ^^oo ^' (3.35) 

This shows that the operator 3" is continuous on the domain D, which is a convex compact set in 
C'^([0,T] X [0, 1])^. It follows then from Schauder fixed-point theorem that 3" has a fixed-point. 
This proves the existence of solutions on the time interval [0, T]. 

3.5. Uniqueness of the solution. Let uq,vq E Lip([0, 1]) be as in ()3.8p - ()3.9p . Assume given 
two pairs {u^,v^), (ti^,w^) E D of solutions of ()3.5p - ()3.7p : that is, if ui and Vr are defined as in 
p.lOp . then li*, i = 1, 2, is a (weak) solution of 

dtu' + \{u\v')dccu' = 0, 
u'{t,0) = ui{t), u'{0,x) = uo{x), 

while f*, i = 1,2, is a (weak) solution of 

dtv"^ + fi{u\v'^)dxv'^ = 0, 
v'{t,l) = Vrit), v'{0,x) = vo{x), 

Let u = u^ -u"^ and v = v^ - v"^. Note that u,v e Lip([0, T] x [0, 1]) = W^'°°{{0, T) x (0, 1)) 
and that u, v fulfill 

dtu + X^d^u + Xd^rU^ = 0, (t,x) E (0,r) X (0, 1), (3.36) 

dtv + fi^d^v + fid^v^ = 0, {t, x) E (0, T) x (0, 1), (3.37) 

u{t,0) = v{t,l) = 0, u{0,x) =v{0,x) = 0, (3.38) 

where A* = A(u*,f*), = fi{u\v'^), and X = X^ — X'^ , fi = ji^ — fj?. 

Multiplying in (|3.36p by 2u, in (|3.37p by 2v, integrating over (0, t) x (0, 1), and adding the 
two equations gives 



{t,x) E (0,T) X (0,1), 



it,x) E (0,r) X (0,1). 



[[""(Olli + 11^(^^)112 + 2 / / {X^udxU + l^^vda:v)dxds + 2 f f (XudxU^ + fivd^v^) dxds = 0. 
Jo Jo Jo Jo 

Using (|3.38p and an integration by parts, we obtain 



t rl 



Jo 

= - f f [{OxX^M"^ + idxH^)\v\'^]dxds+ I [X^\u{s, 1)1'^ - n^\v{s,0)\'^]ds 
Jo Jo 

>- f I [{dxX^m'' + {^x^i^)\v\'']dxds 
Jo Jo 



where we used ()3.3p . On the other hand, since A and /u are M2-Lipschitz continuous on 
[— Ci, Ci\^ , we infer that A* and /u* are 2M2C3-Lipschitz continuous on [0, T] x [0, 1]. In particular, 

WdxX^Woo < 2M2C3, Wdxfl^Woo < 2M2C3 
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and 

|A| < M2{\u\ + \v\), 
lAl < M2{\u\ + \v\). 

This yields 

\2 [ [ {Xud:,u^ + fivd-,v^) dxds\ < 2M2C3 [ [ {\u\ + \v\f dxds . 
Jo Jo Jo Jo 

We conclude that for all t £ (0, T) 

wrnwl + wmwl < 6M2C3 Aiiniii + \\v\\i)ds. 

Jo 

This yields u = v = 0, hy Gronwall's lemma. 

3.6. Finite-time extinction of the maximal solutions. In this section, (u, v) denotes the 
only solution of ()3.5p - ()3.7p in the class D. 

Lemma 3. At time t = T we have 

u{T,x) = v{T,x) = 0, Vxe[0, 1]. (3.39) 

Proof of Lemma 0.' We infer from (j3.6p that 

u(t,0) =v(t,l) =0, \ft>T--. (3.40) 

c 

Thanks to (j3.2p - (|3.3p . we have that 

X{u{t, x),v{t, x)) > c> > -c> n{u{t, x),v{t, x)), V(t, x) G [0, T] x [0, 1]. 

Let (f)'^ (resp. (j)^) denote the flow of X{u,v) (resp. fi{u,v)), and let e'^ (resp. e^) denote the 
corresponding entrance times. (Note that > implies (p^{e^{t,x),t,x) = 1.) Then the 
following holds: 

e''(T,x)>T -- and e^(r,x)>r--, VxG[0,1]. 
c c 

Combining this with (fSlIO]) and ([227]), we obtain ([339]) . □ 
Finally, it is sufficient to extend u and f by for t > T to get a global in time solution. The 
stability property (f3l^ follows at once from as the r.h.s. in (^:2T\) and ([3:22]l 

tend to as (Ci,C2) — >■ (0,0). The proof of Theorem[T]is complete. □ 

4. Finite time stabilization with a control from one side 

In this section, we consider a system of the form 

dtu + X{u,v)d^u = 0, (t,x) G (0,+oo) X (0, 1), (4.1) 

dtv + fi{u,v)d^v = 0, (t,x) G (0,+oo) X (0,1), (4.2) 

u{t, 0) = h{v{t, 0),t), u{0, x) = uo(x), (4.3) 

v{t,l) =Vr{t), viO,x) =vo{x), (4.4) 
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where Vr still solves the ODE 

d 



dt 



Vr{t) 



In (|4.3|) . h denotes some function in Q}{[—Ci,Ci] x 
number Ci > such that, for some time > 0, 



Ksgn{vrit))\vr{t)\\ VriO)=Vo{l). (4.5) 
)nW^'°°{{-C^,C^) X (0,+oo)) for some 



h{0,t)=0 yt>Th. 



We introduce the numbers 

Ci G (0,Cr], 

T :-- 



^1-7 



-+max(Tft,- 
c c (1 — 7)A 



Di 

Ml 
Mo 



max(Ci , 1 1 /i| I L°° ((-Ci ,Ci ) X (o,+oo)) ) ) 

||c^i>^||l°°((-Ci,Ci)x(0,+oo))) 
Pt/l||Loo((_(7i,Ci)x(0,+oo)), 



max 



max 



max( 
max( 



e<^{[-c[,c[]x[-Ci,Ci])A\W\eo{[-c[,c[]x[-c,,c^])), 
\\9ufJ'\\eo([-c[,c[]x[-Ci,Ci])^ \ \9vfJ'\\&\[~c[,c[]x[-Ci,Ci])^ 
\\duM\eo([-c[,c[]x[~Ci,Ci])i \ \9vM\eo{[~-c[,c[]x[-Ci,Ci]))^ 

^ ,C2), 



2TM2 



C2) max(l, Ml) exp(2TM2Q 



3;- 



+ D2. 



Note that, if ||f ||eo([o,T]x[o,i]) ^ C'l, then for all t G (0,T) 

|n(t,0)| < C[ and |aj'u(t,0)| < Dl\^tv{t,^ 
We shall consider the following conditions 

C'^ := max(i(Z)iC^ + i^a), C2) max(l, Mi) exp(2TM2C3) < C3. 

Note that (|4.7p and (|4.8p are satisfied if Ci, C2, and Da are small enough. 
We introduce the set 



(4.6) 



(4.7) 
(4.8) 



D := {{u,v) G Lip([0,r] x [0,1])^ ll^i|leo{[o,T]x[o,i]) ^ ^i^ ll^^lleo{[o,T]x[o,i]) ^ C*!' 

n is Cs-Lipschitz, v is C3-Lipschitz} . 

We pick a pair (no,fo) G Lip([0, 1])^ fulfilling (|3.8p - ()3.9p and the following compatibility 
condition 

uo{0) = h{vo{0),0). (4.9) 
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Let US do some comments about the boundary condition ()4.3p . For a system of conserva- 
tion laws on the interval (0,1), a very general boundary condition at j; = takes the form 
f{u{t,0),v{t,0)) = 0. If duf{uo,vo) 7^ 0, then around {uq,vo) an application of the Implicit 
Function Theorem gives a relation of the form 

u{t,0) = h{v{t,0)) 

with h a smooth function of w in a neighborhood of vq. Assume now that the interval represents 
an edge in a network, and that the left endpoint is a multiple node (i.e. it belongs to at least 
two edges). The contributions of the other edges at this multiple node can be taken into account 
in h through its dependence in t in ()4.3p . 

We are in a position to state the main result of this section. 

Theorem 2. Assume that Ci,C2 and D2 are such that the conditions ()4.7p and ()4.8p are 
satisfied. Then for any pair (uq^vq) G Lip([0, 1])^ fulfilling ()3.8p . ()3.9p and (|4.9p . there exists a 
unique solution {u,v) of (|4.ip - (|4.5p in the class D. Furthermore, the solution is global in time 
with u{t, .) = v{t, .) = for t >T. Finally, if h = h{v), then the equilibrium state (0, 0) is stable 
in Lip([0, 1])^ for ^i^-^M); that is 

ll(^^,'f)llL->{R+;Lip{[0,l])2) as ||(ii0,'f^0)||Lip{[0,l])2 0. (4.10) 

Proof. It is very similar to those of Theorem [TJ If {u,v) G D is given, we define {u,v) = 3'{u,v) 
as follows: u is the weak solution of the system 

jdtU + X{u,v)da;U = 0, \^If^rr^^^n^^ 
\ n^ n\ +\ \ r\ V(t, x) G [0, T] x [0, 1], 

\u{t,0) = n{v[t,0),t), u[0,x) = uq[x), 

and V is the weak solution of the system 

(dtv + n{u,v)d^v = 0, G [0 T] X [0 1] 

]^v{t,l) = Vr{t), v{0,x) = Vo{x), 

Then, using Proposition 12.61 and ()4.7p - (|4.8p . one readily sees that 

ll'"lleo([o,T]x[o,i]) - ^1' ll^llccicTixici]) < Ci, 

u is C3-Lipschitz, hence u is Cs-Lipschitz, (4.11) 
V is Cg-Lipschitz, (4.12) 

so that 3" maps T> into itself. Let us prove that 3" is continuous, D being equipped with the 
topology of the uniform convergence. Consider a sequence {{un,Vn)} C T> and a pair {u,v) G T> 
such that 

max (^\\un - it||eo{[o,T]x[o,i])JI^'n - '^lleo{[o,T]x[o,i])) 0- (4-13) 

Let 

{un,Vn) = 3'iun,Vn) for n > 0, and {u,v) = 3'{u,v). (4.14) 
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We aim to prove that Un ^ u and Vn ^ v uniformly on [0, T] x [0, 1] as n — )• oo. We focus 
on u„, the argument for Vn being the same as those given in Lemma El We consider the same 
(j)n^ (j), e„, e, In, I, Jn, J, Pn and P, as in the proof of Lemma [2j Then 



Un{t,x) = 

and 

u{t, x) 



h{vn{en{t,x),0),enit,x)) if {t,x) e Jn, 

{0,t,x)) if{t,x) £ InUPn 

j h{v{e{t,x),0),e{t,x)) if {t,x) G J, 

\uo{(l){0,t,x)) ii{t,x)£lUP. 



Assume first that {t, x) £ J. Then e{t, x) > and e„(t, x) > for n large enough, by Proposition 
2.21 Since Vn ^ v uniformly on [0, T] x [0, 1] and e„(t,x) — )■ e{t,x), we infer that 



Unit, x) = h{vn{en{t, x), 0), e„(t, x)) h{v{e{t, x), 0), e{t, x)) = u{t, x). 
If now {t,x) G /, one can repeat the argument in Lemma [2] to conclude that 

Un{t,x) = Uo{(pniO,t,x)) ?io((/'(0, t, x)) = u{t,x). 

Thus, Un{t,x) — u{t,x) for (t, x) £ lU J, hence for a.e. (i,x) G [0, T] x [0,1]. We have also 
that Vn{t,x) — )■ v{t,x) for a.e. {t,x) G [0,T] x [0,1]. We infer from the compactness of D 
in e°([0,T] X [0, 1])2 that ('u„,w„) ^ {u,v) in e°([0,r] x [0,1])^. We conclude with Schauder 
fixed-point theorem that 3" has a fixed-point {u,v)gT>, which is a solution of (|4.ip - ()4.5p on 

[o,r]x[o,i]. 

Let us now establish the uniqueness of the solution of ()4.ip - ()4.5p in the class D. Assume given 
two pairs {u^,v^), {u^,v^) G D of solutions of ()4.ip - ()4.5p : that is, with Vr defined as in ()4.5p . f*, 
z = 1, 2, is a (weak) solution of 

^^.^y^'^'^f^'^^rn^S ( ^ (t,x)G(0,T)x(0,l), 

while u*, i = 1, 2, is a (weak) solution of 
dtu' + X{u\v')d,u' = 0, 

MHt,0) = /i(fHt,0),t), u*(0,x) =^0(2;), t ^u,i j X ^u,ij. 

Let u = — and w = — f ^. Note that u,v £ VF^'°°((0, T) x (0, 1)) and that u, v satisfy 

dtu + X^d^u + Xd^u^ = 0, (4.15) 

dtv + ^x^d:,v + fld^v"^ = 0, (4.16) 

u{t,0) = h{v'^{t,0),t) - h{v^{t,0),t), (4.17) 

vit,l)=0, (4.18) 

u{0,x) =v{0,x) = (4.19) 

where A* = X{u\v^), ji^ = /x(u*, u*), and X = X^ — X^ , ji = — jj"^. 
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Multiplying in ()4.15p by 2u and integrating over (0,t) x (0, 1) gives 

||n(t)||2 = _2 [ [ [X^ud^u + Xud^u^]dxds 
Jo Jo 

t rl ft 

[{d,cX^)\u\'^ - 2Xud^u'^]dxds - / X^\u\'^\lds 



^0 ^0 

< 2M2C3 /" [ [\u\'^ + \u\i\u\ + \v\)]dxds 
Jo Jo 



+\\M\&\hCi,ciM-CuC,])Di \v{s,0)\'ds (4.20) 



where we used (|3.3|) . 

Multiplying in ()4.16p by 2v and integrating over (0, t) x (0, 1) gives 

~ ~^ [fi^vdxV + fivdxV^]dxds 

Jo Jo 

= [ f [{dxiJ^)\vf -2fivdxv^]dxds + [ ij^\v{s,0)fds 
Jo Jo JO 



< 2M2C3 [ [ [\v\^ + \v\{\u\ + \v\)]dxds - c [ \v{s,0)\''^ds (4.21 
Jo Jo Jo 



where we used p.3p again. Let us introduce the energy 

E{t) = \\u{t)\\^ + \\X\\eoi[^c[,c[M^CuC^])^Mt)f- 
Combining (fOOj) with ([OT]) yields 

E{t) <C I E{s)ds, 
Jo 

for some C depending only on T), so that -E = 0, by Gronwall's lemma. This proves the 
uniqueness. For the extinction time, we notice that from the proof of Theorem [1] 

1 (7I-7 

v(t, x) = 0, for - + - — <t<T, < X < 1. 
c (1 — 7)a 

Combined with (|4.6p . this yields 

/ ^ (ji-j \ 
uit,0) = h{vit,0),t) = 0, for max ( ^fe, - + ^^^^ j <t< T. 

Using (jS.Sp . we conclude that 

u{T,x) = VxG[0, 1]. 

Assume now that h = h{v), i.e. D2 = 0. The stability property ()4.10p follows at once from ()4.6p 
and (1121)- (HSU, as C[ < max(l,Di)Ci and iC%,C^) (0,0) as (Ci,C2) (0,0). The proof of 
Theorem [2] is complete. □ 
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5. Application to the regulation of water flow in channels 



In this section, we investigate the regulation of water flow in a network of open horizontal 
channels. We assume that the channels have a rectangular cross section and that the friction on 
the walls can be neglected. In this context, the flow of the fluid can be described in a satisfactory 
way by the shallow water equations (also called Saint- Venant equations) (see l20j). The control 
in feedback form is applied at the vertices of the network, which is assumed to be a tree. 

We introduce some notations needed in what follows (we follow closely [S]). Let T be a tree, 
whose vertices (or nodes) are numbered by the index n £ = {1, A^}, and whose edges are 
numbered by the index i £ 3 = {!,...,/} with I = N — 1. We choose a simple vertex, called 
the root of T and denoted by 3J, and which corresponds to the index n = N. We choose an 
orientation of the edges in the tree such that 3? is the "last" encountered vertex. It is similar to 
those of a fluvial network in which each edge stands for a river, and 3? indicates the place where 
the last river enters into the ocean. 

We denote by li the length of the edge with index i. Once the orientation is chosen, each 
point of the i-th edge is identified with a real number x £ [0, k]. The points x = and x = k 
are termed the initial point and the final point of the i-edge, respectively. 

Renumbering the edges if needed, we may assume that the edge with index i has as initial 
point the vertex with the (same) index n = i for all i £ 3. 

We denote by J„ C J, n = 1, A, the set of indices of those edges having the vertex of index 
n as one of their ends. Let 



Note that Ei^i = for all i £ 3, and that eN-i,N = 1. A node with index n is said to be simple 
(resp. multiple) if i^{3n) = 1 (resp. #{3n) > 2). The sets of indices of simple and multiple nodes 
are denoted by J^s and J^m, respectively. The depth of the tree is the greater number of edges 
in a path from one simple node to 31. 





1 



if the vertex with index n is the initial point of the edge with index i; 
if the vertex with index n is the final point of the edge with index i. 




13 



umn 



Figure 2. A tree with 14 nodes, a depth equal to 5, with simple nodes 'Ns = 
{1, 2, 4, 5, 6, 7, 8, 14} and multiple nodes T^m = {3, 9, 10, 11, 12, 13}. 



FINITE- TIME STABILIZATION OF SYSTEMS OF CONSERVATION LAWS ON NETWORKS 23 

Pick any channel represented by (say) the i-th edge of the tree, which is identified with the 
segment [0,/j]. Then the shallow water equations read 

dtHi + d^iHiVi) = 0, t>0, 0<x<li, (5.1) 

dtVi + d^{^ +gHi) = 0, t > 0, < X < k, (5.2) 

where Hi{t, x) (resp. Vi{t, x)) is the water depth (resp. the water velocity) along the i-th channel, 
and g is the gravitation constant. The equations ()5.ip - (|5.2p have to be supplemented with some 
initial conditions 

Hi{0,x) = Hi^o{x), mO,x) = Vi,o{x), 0<x<li (5.3) 

and with two boundary conditions. In general, there are at the two ends of the channel (i.e. at 
X = and at x = /j) some hydraulic devices to assign the values of the flow rate. Recall that 
the flow rate is defined along the channel as 

Qi{t,x) := Hi{t,x)Viit,x). 

At any multiple node n S J^m, the equation of conservation of the flow 

^(-l)^-"Q,(t,ei,„Z,) = (5.4) 

has to be taken into consideration. It yields a boundary condition (coming from the physics) 
in which no control applies. Let Iq £ 3n be the only index such that ejo,n = 0, namely io = n. 
Then (|5.4p can be written 

Qio(*>o)= Qiit^^i)- (5-5) 

Thus, the flow rate may be controlled at the final points of the edges of indices i ^ Iq, while it 
is prescribed by ()5.5p at the initial point of the edge of index io. 

We aim to stabilize the system around some equilibrium state, represented by a sequence 
{{H* , V*)}i<i<i of pairs of positive numbers. Let Q* = H*V* . For ()5.4p to be valid as t — )• oo, 
we impose that 

J^(-1)^-"Q* = 0, VnGTlM. (5.6) 
Introduce the characteristic velocities 



Vi - ^/gHi, (5.7) 

\i = y^ + v^ (5.8) 

and the Riemann invariants (see \13\ [20] ) 

Ui = Vi + 2y^,-{V* + 2y^), (5.9) 

Vi = Vi.-2y^i-{V* -2y^). (5.10) 

We shall assume thereafter that the flow is subcritical or fluvial; that is, the characteristic 
velocities are of opposite sign 

/ij < < Aj. 
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Clearly, this holds if 

o<v*< (5.11) 

and max{\Hi — H*\, |Vi — V^*|) is small enough. From now on, we assume that (j5.1ip holds for 
all i G 3, and we pick a number c > such that 



vW-^i*>2c, Vie J. (5.12) 
Note that ()5.9p - ()5.10p may be inverted as 



1 



Hi = ^^H* + ^{u,-v,)j , (5.13) 

Vi = V* + ^{ui + Vi). (5.14) 
Substituting the values oi Hi,Vi in ()5.9p - (|5.10p yields 

fii = V* -^/^+^{ui + 3v^), (5.15) 

Xi = V* + y^+^i3u,+v^). (5.16) 

Combined with (j5.12p . this shows that 

max(|iij|, |uj|) < c =^ fj.i < —c < c < Aj. 

The shallow water equations ()5.ip - ()5.2p . when expressed in terms of the Riemann invariants Ui 
and Vi, read 

dtUi + Xi{ui,Vi)dxUi = 0, t > 0, < X < li, (5.17) 

dm + iJLi{ui,Vi)dxVi = 0, t > 0, < X < /j. (5.18) 

Let us now turn our attention to the boundary conditions. Consider first a boundary condition 
associated with an active control, e.g. 

^^Vi{t,li) = -Ksgn{vi{t,li))\v^{t,li)\^ ■ (5.19) 

In practice, one would like to assign the value of Qi{t, li) = Hi(t, li)Vi{t, li) by using the output 
Hi{t,li) only. Using (j5.1(jp . it is sufficient to set 



Q,{t, k) = Hi{t, li) [vi{t, k) + 2^/gH,it,li) + - 2^^] , (5.20) 
where Vi solves ()5.19p together with the initial condition 

Vi{0, k) = Vi{Q, li) - 2^gHi{0,li) - V* + 2y^. (5.21) 
For a control applied to the initial point of the i-edge, we set 

Qi{t, 0) = H,{t, 0) {ui{t, 0) - 2^gHi{t,0) + V* + 2 yW) > (5-22) 
where Ui{.,0) solves 

j^Ui{t,0) = -Ksgn{ui{t,0))\ui{t,0)\\ (5.23) 



tx,(0,0) = V-(0,0) + 2^5^^,(0,0) - V* - 2y^. (5.24) 
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Consider next a boundary condition without any active control. For a simple node n G Jsf^ and 
the corresponding edge z E J„, if Si^n = (i.e. the node n is the initial point of the edge i), then 
n = i and a natural boundary condition at the node n is given by the relation 

Qiit,0) = Q*, (5.25) 



that is 
where 

Since 



Fi{u^{t,0),vi{t,0))=0 (5.26) 

a TP -1 T^* 

F,(0,0) = and ^(0,0) = -/^(/^+ -^) > 

it follows from the Implicit Function Theorem that there exist a number 6i > and a function 
hi G C^(M) with /ij(0) = such that for max(|?x|, |z;|) < 6i, 

Fi{u,v) = <;=^ u = hi{v). 

Thus (|5.26p may be written, at least locally, in the form 

u,{t,{)) = hi{vi{t,0)). 

Finally, for a multiple node n G ?^m, if io £ 3n is the only index such that ejg^ri = (i.e. io = n), 
then ()5.5p may be written 

Fi„(f/,„(t,0),^;io(t,0),[/(t),y(t)) = 
where ?7(t) = (uj(t, /i))jea„,i^jo, V{t) = {vi{t,li))i(zj^^i^i^ and 

Note that, by ([SSI), -FjoCO' 0' 0' 0) = and 

BF 1 I I V* 

'■"-(0,0,0,0) = -,/^(,/^+^)>0. 



dui,' ' ' ' ' 2V *o 

We may pick a number > and a function i/jg of class around such that, if \ui\ < 
and It'll < (5io for all i G J^, we have 

Fio {uio ,Vig,U,V) = <^ Ui,) = (wio ,U,V). 

Replacing {ui,Vi) by {ui{t,li),v{t,li)) for i / io in C/, F, we see that ()5.5p may be written, at 
least locally, in the form 

u,,{t,0) = hi,{v,,{t,0),t) (5.27) 

where /lig G e^(M2) and /ij(,(0,t) = if Ui{t,li) = Vi{t,li) = for ah i G J„ \ {io}- 

We are in a position to state our results for the regulation of water flow in channels. Consider 
first one channel (N = {1, 2}, J = {!}) represented by the segment [0, ^i]. 



26 VINCENT PERROLLAZ AND LIONEL ROSIER 

Theorem 3. (Two boundary controls) Assume that (jS.lip holds for i = 1, and pick any c > 
as in ()5.12|) . Then there exists a number 6 > such that for all {Hifi,Vifi) G Lip([0, Zi])^ with 

maxdli^i^o - -f^i lli4^i.°°(o,Zi)) ll^i,o - ^1*11^1. °°(o,h)) < (5.28) 

there exists for any T > a unique solution {Hi,Vi) G Lip([0, T] x [0,/i])^ of (j5.ip - (|5.3p and 
(|5.19p - (j5.24p . Furthermore, there exists a function t*{H^,Vi,c, 5, Kjj) with lim^^Q t* = sttc/i 
that 

Hi(t,x) = m, Vi(t,x) = V{ t>- + t*,x£(0,li). (5.29) 

c 

Finally, the equilibrium point {H^,Vi) is stable in Lip([0, 1])^ for the system (|5.ip - (|5.3p and 

Proof. Noticing that the map : {Hi, Vi) — )• {ui,vi) defined along ()5.9p - ()5.10p is locally around 
(i?i*, Vi*) a diffeomorphism of class 6°°, the condition ^^2Eh implies ([331)- ([221) for Ci and C2 as 
in Theorem [1] (applied actually on the interval (0, ^i) rather than (0, 1)), provided that 5 < (^0 is 
small enough. We modify the functions ^i{u,v) and Xi{u,v) outside [— c, c]^ so that 

fii{u,v) < —c < c < Xi{u,v), {u,v) G M^. 

Let {ui,vi) be the solution given by Theorem[Tl and let {Hi, Vi) := Q~^{ui,vi). If Ci is chosen 
sufficiently small, then we infer from ()3.2ip that 

max(|ui(t, a;)|, |t;i(t, x)|) < Ci < c, t > 0, < x < /i, 

mscK{\Hi{t,x) - Hl\,\Vi{t,x) -Vi\) < 5, t>0, 0<x</i. 

It follows that for all T > 0, {Hi, Vi) G Lip([0, T] x [0, li\f is a solution of ([SID-dLSl) and (I5J9D - 
([Qi]) such that (fog]) holds with = Cl^'^ /{{I - ^)K). Note that the range of Ci in Theorem 
[1] depends on H^, V{ and c through the constants Mi and M2, and that (Ci,C2) — )• (0,0) as 
(5 — )■ 0. Thus t* — 7- as (5 — )• with and 7 kept constant. The uniqueness of {Hi, Vi) in the 
class Lip([0,r] X [0,li])^ for all T > follows at once from those of {ui,vi) in the same class, 
as stated in Theorem [TJ The stability property follows from ()3.18p . □ 

If the control is active at one endpoint of the channel only, a finite-time stabilization may be 
derived as well. 

Theorem 4. (One boundary control) Assume that (jS.lip holds for i = 1, and pick any c > 
as in (j5.12p . Then there exists a number 6 > such that for all {Hifi,Vi,o) G Lip([0, Zi])^ with 

max(||iJi,o - -ffillH/i.°°(0,«i) Jl^i,0 - V'i*||vyi,oo(o,i,)) < S, (5.30) 

^1,0(^1)^1,0(^1) =^i*^^r, (5.31) 

there exists for anyT >0 a unique solution {Hi,Vi) G Lip([0, T] x [0, h])^ of ([SJD-dSS]), ([5T9]1 - 
(|5.2ip and (|5.25p . Furthermore, there exists a function t*{H*,Vi,c,6,K,j) with lim5_>.o = 
such that 

11 

Hi{t,x) = Hl, Vi{t,x) = Vi t > — + r, X G (0,/i). (5.32) 

c 

Finally, the equilibrium point {HI, Vi) is stable in Lip([0, 1])^ for the system ([539]) - 
(fOT]) and ([5^25]) . 
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Proof. It is sufficient to proceed as for the proof of Theorem [3l and to use Theorem [2] (on the 
domain [0, li] and with a control active at the final point only). □ 

A direct application of Theorem S] gives the following result for a chain of two channels 
(?sf = {1, 2, 3}, J = {1, 2}), for which there is no active control at the internal node. 

Corollary 1. (Two channels and two controls) Assume that (|5.1ip holds for i = 1,2 with 
Qi = Qh ^''^d pick any c > as in (j5.12p . Then there exists a number 5 > such that for all 
iHi,o,Vi,o,H2,o,V2,o) e Lip([0,Zi])2 X Lip([0,Z2])' with 

max{\\Hifi - 1,4/1,00(0,;^), 11^^,0 - l^i*||iyi.°°(o,«,)) < ^ = 1,2, (5.33) 

i^i,o(^i)V^i,o(/i) = Ql = Q*2 = ^2,o(0)F2,o(0), (5.34) 

there exists for any T > a unique solution {Hi,Vi) G Lip([0,T] x [0,/i])^ of (|5.ip - (|5.3p for 
i = l,2, <^M)-§2Ii) for i = 2, ^^-^M) for i = I, and 

Qi{tM) = Q*i = Ql = Q2{m, t>0. 

Furthermore, there exists a function t*{H^, Vj^*, V^*, c, 6, K, 7) with lim^^^o t* = such that 

Hi{t,x) = H*, Vi{t,x) = V* t> ^™^^(^i'^2) ^^.^ xe{0,li), i = l,2. (5.35) 

c 

We can extend the above results to a network of open channels which is a tree. We assume that 
the incoming flows can be controlled at each multiple node (the outgoing flow being uncontrolled 
and deduced from the conservation of the flows). In terms of Riemann invariants, for the edge 
with index i, the function Vi is controlled at x = li according to ()5.19p . while the function Ui 
is controlled at j; = according to ()5.23p only if the initial point of the edge is a simple node 
(otherwise, Ui{t,0) is given by ()5.27p ). 

The main result of this section is the following 

Theorem 5. (Network of open channels) Consider a tree with N nodes and I = N — 1 edges. 
Assume that (|5.1ip holds for i = that (|5.6p holds, and pick any c > as in (|5.12p . 

Then there exists a number 6 > such that for all (-ffi,o, ^1,0; •••) Hj^, V/,o) £ Lip([0) ^1])^ x • • • x 
Lip([0,//])2 with 

max{\\Hifi - -ffriliyi.->(0,«0' H^^-O - V*\\w^.^(o,h)) < ^ = 1, --Ji (5-36) 

^^n,o(0,0)K,o(0,0) = Yl H,,oiO,h)Vi,oiO,h), VnGNM, (5.37) 

there exists for any T > a unique function {Hi,Vi, Hf,Vi) S Lip([0,T] x [0,/i])^ x ••• x 
Lip([0,r] X [0,li]f such that, for all i = 1,...,/, (I5l])-(I53D and (l5T9]l - (i5:2T|) hold, and (15:22]) - 
()5.24p hold if the initial point of the i-th edge is simple, while (|5.5p holds if the initial point of 
the i-th edge is multiple. Furthermore, there exists a function t*{Hl,V{ , Hj ,Vj ,6,c, K,j) 
with lim5_j.o t* = such that 

H,it,x) = H*, V^{t,x) = V* Pmaxi<,<,/, x e {0,li), i = !,...,!, (5.38) 

where p denotes the depth of the tree. Finally, the equilibrium state {H* ,V*)i<i<i is stable in 
Lip([0, /i])^ X • • • X Lip([0, //])^ for the system. 
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Proof. The proof is done by induction on the number of edges / > 1. For 1 = 1, the result was 
already proved in Theorem[3l Note that the norm ||(^ii,fi)||Loo(R+;Lip([o,«i])2) in Theorem[3]is as 
small as desired if 6 is small enough. Let I > 2, and assume the result true for any tree with at 
most I — 1 edges, with the norms IK^^i, i'i)llL°°(R+;Lip([o,/i])2) the edges of the tree as small as 
desired if 6 is small enough. Pick any tree with I edges. Recall that the root 31 is the node with 
index A^, and that it is the final point of the edge of index I = N — 1. Denote by 3?' the initial 
point of the edge of index /, i.e. the node of index N — 1. Let k = #{3n-i), and let us denote 
by Ti, Tfc_i the subtrees of T with 31' as root. (3? does not belong to any of them.) Note 
that the subsystem associated with any subtree Tj is decoupled from the other subtrees and 
from the last edge of index /. An application of the induction hypothesis on each subtree Tj, 
1 < i < k — 1, yields the existence (and uniqueness) of the functions (Hi, Vi) for i = 1, I — 1. 
Next, the existence and uniqueness of {Hj,Vi) follows at once from Theorem [21 Indeed, the 
constant D2 in Theorem [5] may be taken as small as we want if 6 is sufficiently small, for the 
quantities \\dtUi{., li)\\oo and \\dtVi{., li)\\oo for i G Jat-i \ {N — 1} may be taken arbitrarily small 
by the induction assumption. Furthermore, the norm ||(u/,i^7)||Loo(K+.Lip([o^;^])2) tends to with 
6, by ()4.1ip - ()4.12p . The condition ()5.38p is obtained by an obvious induction on the depth of 
the tree. □ 



Appendix: Proof of Proposition 12.11 

First, we introduce some extension operator n which maps a function a : [0, T] x [0, 1] 
to a function a = 11(a) :M? ^M. defined as follows: 



for < t < T 

a(t, x) 



a{t,x) if X e [0, 1], 

a{t,2-x) if X e [1,2], 



and a{t,x) is 2-periodic in x (i.e. d{t,x + 2) = d(t,x)); 

• for t > r, d{t, x) = d{T, x) for all x G M; 

• for t < 0, d{t, x) = a(0, x) for all x G M; 

It is easy to see that 11 is a (linear) operator from C''([0,T] x [0,1]) to C'^(M^) (resp. from 
L°°(0,r;Lip([0, 1])) to L~(M;Lip(M)) such that 

l|n(a)||Loo(]R2) = ||a||eo([o,T]x[o,i])i (5.39) 

l|n(a)||ioo(iR.Lip(]8)) = ||a||Loo(0,T;Lip([0,l]))) (5.40) 

n(a)(t, x) > c V(t, x) G if a{t, x) > c V(t, x) G [0, T] x [0, 1]. (5.41) 

Let a fulfill (j2.2p . and let (p (resp. (j)) denote the flow associated with a (resp. with d = vr(a)). 
Then 6 is defined in M.^, and 



(p{s, t, x) = (/)(s, t, x) V(s, t, x) G Dom (j). 



(5.42) 
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Thus it is sufficient to prove that cj) is X-Lipschitz on [0, T]^ x [0,1]. To this end, pick any 
{si,tuxi),{s2,t2,x2) £ [0,r]2 X [0,1]. Then 

\(t>{si,ti,Xi) - 4>{s2,t2,X2)\ 

< \4>{si,ti,Xi) - (f>{s2,ti,Xi)\ + \4>{s2,ti,Xi) - 4>{s2,t2,Xi)\ + \4>{s2,t2, Xi) - (f>{s2,t2, X2)\ 

First, 

2 _ l'S2 

ds4>{T,ti,xi)dT\ = \ / a(r, <^(t, ti, xi))dr| < I |a||eo([o,T]x [0,1]) I Si - •S2I, (5.43) 

Sl J Si 

where we used (|5.39p . For I2, we notice that for all s 

\ds[4>is,ti,xi) - 4){s,t2,xi)]\ = \a{s,4>{s,ti,xi)) - d{s,(f>{s,t2,xi))\ 

< L\(f>{s,ti,xi) - (f>{s,t2,Xl)\ 

where we used ()5.40p . Gronwall's lemma combined to the estimate for Ii yields then 

\4>{s2,tl,Xi) - 4>{s2,t2,Xi)\ < \<p{t2,ti,Xi) - 4>{t2,t2,Xi)\e^^'^-^^^ 

< \4>(.t2,ti,xi)-4>{ti,ti,xi)\e^'^ 

< l|a||eo{[o,r]x[o,i])«^'^l*i - *2|- (5.44) 

Finally, for I3, we notice that for all s 

\ds[4>is,t2,Xi) - 4>{s,t2,X2)]\ = \d{s,4>{s,t2,Xi)) - d{s,(f>{s,t2,X2))\ 

< L\4>{s,t2,Xi) - (f){s,t2,X2)\, 

which, combined with Gronwall lemma, yields 

\4>is2,t2,Xi)-^{s2,t2,X2)\ < e^^'^-^^^\xi - X2\ < e^^\xi - X2\. (5.45) 

Then (USD follows at once from (fSlIHD - ffOSD . 
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